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ABSTRACT 

We present the selection and classification of over a thousand ultraviolet (UV) variable sources 
discovered in ^ 40 deg^ of GALEX Time Domain Survey (TDS) NUV images observed with a 
cadence of 2 days and a baseline of observations of 3 years. The GALEX TDS fields were designed 
to be in spatial and temporal coordination with the Pan-STARRSl Medium Deep Survey, which 
provides deep optical imaging and simultaneous optical transient detections via image differencing. 
We characterize the GALEX photometric errors empirically as a function of mean magnitude, and 
select sources that vary at the Su level in at least one epoch. We measure the statistical properties 
of the UV variability, including the structure function on timescales of days and years. We report 
classifications for the GALEX TDS sample using a combination of optical host colors and morphology, 
UV light curve characteristics, and matches to archival X-ray, and spectroscopy catalogs. We classify 
62% of the sources as active galaxies (358 quasars and 305 active galactic nuclei), and 10% as variable 
stars (including 37 RR Lyrae, 53 M dwarf flare stars, and 2 cataclysmic variables). We detect a large- 
amplitude tail in the UV variability distribution for M-dwarf flare stars and RR Lyrae, reaching up to 
I Am| — 4.6 mag and 2.9 mag, respectively. The mean amplitude of the structure function for quasars 
on year timescales is 5 times larger than observed at optical wavelengths. The remaining unclassified 
sources include UV-bright extragalactic transients, two of which have been spectroscopically confirmed 
to be a young core-collapse supernova and a flare from the tidal disruption of a star by dormant 
supermassive black hole. We calculate a surface density for variable sources in the UV with NUV < 23 
mag and \Am\ > 0.2 mag of ~ 8.0, 7.7, and 1.8 deg~^ for quasars, AGNs, and RR Lyrae stars, 
respectively. We also calculate a surface density rate in the UV for transient sources, using the 
effective survey time at the cadence appropriate to each class, of ~ 15 and 52 deg~^ yr~^ for M 
dwarfs and extragalactic transients, respectively 
Subject headings: ultraviolet: general — surveys 



1. INTRODUCTION 

Unlike the optical. X-ray, and 7-ray sky, which have 
been systematically studied in the time domain in the 
search for supernovae (SNe) and gamma-ray bursts 
■ (GRBs), the wide- field UV time domain is a relatively 
unexplored parameter space. The launch of the GALEX 
satellite with its 1.25 deg diameter field of view, and lim- 
[ iting sensitivity per 1.5 ks visit of ^23 mag in the FUV 
\ (Ac ff = 1 539 A) and N UV (Acff = 2316 A) (Martin_eLaD 
[2OO5I: iMorrissev et al.l2007D , enabled the discovery of UV 
variable sources in repeated observations over hundreds 
of square degrees for the first time. 
The UV waveband is particularly sensitive to hot 
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(w 10* K) thermal emission from such transient and 
variable phenomena as young core-collapse SNe, the in- 
ner regions of the accretion fiow around accreting su- 
permassive black holes (SMBHs), and the flaring states 
of variable stars. The characterstic timescales of vari- 
ables and transients in the UV range from minutes to 
years. M dwarf flare stars have strong magnetic activ- 
ity that manifests itself in fia res of thermal UV em ission 
on the timescale of minutes ()Kowalski et al.l [20091) . RR 
Lyrae stars have periodic pulsations which drive tem- 
perature fiuctuations from ~ 6000 to 8000 K that cause 
periodic variability in the UV on a timescale of ^ 0.5 d 
(jWheatley et al.l[2012l ). Core-collapse supernovae (SNe) 
remain bright in the UV for hours up to several days 
following shock bre akout, dep ending on the radius of the 



proge nitor star ( N akar fc Sari 2010; Rabinak fc Waxman 
'2011') , and the presence of a dense wind (lOfek et al 



2010t iChevalier fc IrwinllMlL iSvirski et al. l l2012D . Ac 
tive galactic nuclei (AGN) and quasars demonstrate 
stronger variability with decreasing wavelength and 
longer timescales of years (jVanden Berk et al.l[200l . 

UV variability studies of GALEX data observed as 
part of the All-Sky, Medium, and Deep Imaging base- 
line mission surveys (AIS, MIS, DIS) from 2 003 to 2007, 
yield ed the detection of M-dwarf fiare stars ("Welsh et al " 



2007 ) , RR Lyrae stars, AGN, and quas ars (Welsh et al 



2005HWheatlev et al.ll2008H Welsh et alllMllI) . and flares 
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from the tidal disrupti on of stars around dormant su - 
permassive black holes (jCezari et al.l[2(M [2008al 120091 ) . 
Serendipitous overlap of 4 GALEX DIS fields with the 
optical CFHT Supernova Legacy Survey, enabled the ex- 
traction of simultaneous optical light curves from image 
differen cing for 2 of the tidal disruption event (TDE) can- 
didates ([Gezari et al.l2008al) , and enabled the association 
of transient UV emission with two Ty pe IIP supernovae 
(SNe ) within hours of sho ck breakout ([Schawinski et all 
[200l iGezari et al.] |2008bl. Chance overlap of GALEX 
observations with the survey area of the optical Palomar 
Transient Factory (PTE) detected a Type Iln SN whose 
rising UV emission over a few days was interpreted as 
a delayed shock breakout through a dense circumstellar 
medium ()Ofek et al.ll2010l ). 

Motivated by the promising results from the analy- 
sis of random repeated GALEX observations, and the 
demonstrated value of overlap with optical time domain 
surveys, we initiated a dedicated GALEX Time Domain 
Survey (TDS) to systematically study UV variability 
on timescales of days to years with multiple epochs of 
NUV images observed with a regular cadence of 2 days. 
The GALEX TDS fields were selected to overlap with 
the Pan-STARRSl Medium Deep Survey (PSl MDS) 
(jKaiser et al.|[2010l ). GALEX TDS and PSl MDS are 
well-matched in field of view, sensitivity, and cadence 
(shown in Table [T]) . Here we present the analysis of 42 
GALEX TDS fields which intersect with the PSl MDS 
footprint, for a total area on the sky of 39.91 deg^, which 
were monitored over a baseline of 3.32 years (February 
2008 - June 2011). In this paper we use PSl MDS 
deep stack catalogs to characterize the optical hosts of 
GALEX TDS sources. Simultaneous UV and optical 
variability of GALEX TDS so urces culled from matches 
with the PSl transient alerts (|Huber et al.ll201ll) will be 
presented in future papers. 

The paper is structured as follows. In [J2] we describe 
the GALEX TDS survey design, and in 331 we describe 
our statistical methods for selecting UV- variable sources 
and characterizing their UV variability. In §4] we de- 
scribe the multiwavelength catalog data used to identify 
the hosts of the UV-variable sources, including archival 
optical data, a deep PSl MDS catalog, and archival red- 
shift and X-ray catalogs. In ij5]we describe our sequence 
of steps for classifying the GALEX TDS sources, in [JS]we 
summarize our classification results, and the UV variabil- 
ity properties of our classified sources. In ij7]we conclude 
with implications for future surveys. 

2. GALEX TDS OBSERVATIONS 

GALEX TDS monitored 6 out of 10 total PSl 
MDS fields, with 7 GALEX TDS pointings (labeled 
PS-fieldnameJAOS pointing) at a time to cover the PSl 
7 deg^ field of view. During the window of observing 
visibility of each GALEX TDS field (from 2-4 weeks, 
1 — 2 times per year) , they were observed with a cadence 
of 2 days, and a typical exposure time per epoch of 1.5 
ks (or a 5a point-source limit of uiab ~ 23.3 mag), with 
a range from 1.0 to 1.7 ks. The NUV detector devel- 
oped a problem on 2010 May 4 during observations of 
PS_ELAISN1, and so we do not include epochs observed 
between this time and when the the instrument was fixed 
on 2010 June 23 in our analysis. Figure [T] shows the posi- 
tion of the GALEX TDS fields relative to the PSl MDS 
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Fig. 1. — GALEX TDS 1.1 deg diameter field pointings shown in 
blue, and the PSl MDS 3.5 deg diameter field pointings shown in 
green. Orange hatched regions show the coverage of SDSS in the 
optical, red hatched regions shows the coverage of SWIRE in the 
optical, and cyan rectangles indicate the coverage of the CFHTLS 
Deep (solid lines) and Wide (dashed lines) surveys in the optical. 
Light grey regions show the coverage of XMM-Newton X-ray ob- 
servations, and dark grey regions show the coverage of Chandra 
X-ray observations. 

fields, and Table [5] lists the R.A. and Dec of their cen- 
ters, the Galactic extinction {E(B — V)) for each field 
from the lSchlegel et al.l (|1998[ ) dust maps, and the num- 
ber of epochs per field. The median number of epochs 
per field is 24. PS_CDFS_MOS00 is an exception with 
114 epochs, because it was monitored with a rapid ca- 
dence {\At\ ~ 3 hours) over a period of 10 days in 2010 
November. Some offsets in the GALEX pointings from 
the footprint of the PSl MDS fields were necessary in 
order to avoid UV-bright stars in the field of view that 
would violate the detector's bright-source count limits. 

Figure [2] shows the temporal sampling of GALEX TDS 
observations in the NUV in comparison to the PSl MDS 
observations in the gpi, rpi, ipi, zpi, and ypi bands 
from February 2008 to June 2011. PSl began taking 
commissioning data of the MDS fields in May 2009, but 
did not begin full survey operations until a year later. 
The GALEX FUV detector became non-operational in 
May 2009, and so we only include NUV images in our 
study. 

3. STATISTICAL MEASUREMENTS 

3.1. Selection of Variable Sources 

Since most galaxies are unresolved by the GALEX 
NUV 5.3 arcsec full-width at half-maximum (FWHM) 
point spread function (PSF), we can use simple aper- 
ture photometry instead of image differencing to mea- 
sure variability. We create a master list of unique 
source positions froni the pipeline-generated catalogs 
(jMorrissev et al.l [20071 ) for all the individual epochs, as 
well as deep stacks of all the epochs, using a cluster- 
ing radius of 5 arcsec. This radius is chosen such that 
for the typical astrometry error of GALEX of 0'.'5, the 
Bayesian probability that the match is real is larger 
than the Bayesian probab ility that the match is spurious 
(|Budavari fc Szaiavll2008[ ) . The final master list includes 
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Fig. 2.— Dates of GALEX TDS NUV observations compared 
to PSl MDS observations in the ^pi, rpi, 2pi, 2pi, and ypi bands. 
Dotted lines show yearly intervals. 

419,152 sources. 

In order to select intrinsically variable sources in our 
survey, we first need to characterize the photomet- 
ric errors. Although the GALEX images are Poisson- 
limited, the Poisson error underestimates the total error 
in the GALEX catalo g magnitudes by a factor of ^ 2 
(jTrammel l et al.ll20Q7D . This discrepancy is attributed 
to systematic errors such as uncertainties in the detector 
background and flat-field. Thus, we measure the photo- 
metric error empirically by calculating the standard devi- 
ation of aperture magnitudes in bins of mean magnitude, 
(m). We only include objects in the pipeline-generated 
catalogs that are detected in all or >10 epochs. In each 
bin of objects with {m)i — (^^^if^i,k)/'^ (each bin 
typically has A'' = 50 to 1000 sources), we calculate for 
each epoch of a total of n epochs. 



a{{m)^ k) 



\ 



N 



N -1^ 

1=1 



{rrii 



(1) 



where k is the magnitude (in the AB system) given 
by m,,fe = -2.51og(/6) + zp + C^^ -^.2E(B - V), fe 
is the background-subtracted fiux in a 6 arcsec radius 
aperture, zp = 20.08, the apertur e correction is Cap = 
—0.23 mag (jMorrissev et all [20071 ), and we correct for 
Galactic extinction using the values for E{B — V) listed 
in Table [21 We use 3cr clipping to remove outliers in the 
calculation of a{{m), k) which can arise from artifacts. 

The astrometric precision depends on the signal-to- 
noise of the source, thus we also empirically measure a 
magnitude-dependent clustering radius. We do so by 
measuring the cumulative distribution of spatial sep- 
arations between the position in each epoch and the 
mean position for sources in bins of (m), and record 
d95{{m),k), the value for which 95% of sources have a 
separation less than or equal to that value. The re- 
sulting value for dgs is a strong function magnitude, 
increasing from ^ 1 arcsec for (m) = 18.0 mag to ^ 
4 arcsec for (m) — 23.0 mag. Figures |3] and 0] show 
(j{{m),k) and dg5({m),k) for an example GALEX TDS 




Fig. 3. — Empirical determination of la photometric errors 
as a function of mean magnitude from the standard deviation 
of sources detected by the pipeline catalogs in > 10 epochs for 
one of the GALEX TDS fields PS.COSMOS.MOS23. Solid green 
line shows a quadratic fit to the error function for the field 
PS_COSMOS_MOS23, and dashed blue line shows a quadratic fit 
to the median error function for all of the GALEX TDS fields. 
Solid red line shows the expected la Poisson error, which underes- 
timates the total photometric error by a factor of ~ 2. 




Fig. 4. — Maximum spatial separation from the mean for 95% of 
the sources as a function of mean magnitude from the cumulative 
distribution of sources detected by the pipeline catalogs in > 10 
epochs for one of the GALEX TDS fields PS.COSMOS.MOS23. 
Solid green line shows a quadratic fit to the error function for the 
field PS.COSMOS.MOS23, and dashed blue line shows a quadratic 
fit to the median distance function for all of the GALEX TDS 
fields. Note that due to systematic differences in the PSF be- 
tween fields, the distance function radius for PS.COSMOS.MOS23 
is up to ~ 1 arcsec smaller than the median for all fields, since for 
(m) < 20 mag, dgs ~ 2 arcsec in the PSJCMMLSS, PS.CDFS, 
PS.GROTH, and PS.ELAISN1 fields and dgs ~ 3 arcsec in the 
PS.VVDS22H fields. 

field PS.COSMOS.MOS23, a quadratic fit to the me- 
dian function for all epochs in that field, and the median 
function fit over all fields. 

In our master list of source positions, we include all 
sources detected, including sources detected in only one 
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epoch, and fix the centroid to the epoch for which 
the source is detected with maximum flux. We mea- 
sure forced aperture magnitudes at the positions of each 
source in epochs where the source was not detected by 
the pipeline or the spatial separation of the matched 
source is greater than (i95((m},fc). When the aper- 
ture magnitude is fainter than miim in an epoch, it is 
flagged as an upper limit and replaced with mum, where 
miim = -2.5 log(5^/ (BskyiVpix/roxp) + zp + Cap, where 
BsVy = 3 X 10~^ counts s~^ pixel"^, iVpix = IGtt and Toxp 
is the exposure time of that epoch in seconds. 

We select sources that have at least one epoch for which 
\mk — (to) I > 5a{{m),k), where (m) is calculated only 
from epochs that have a magnitude above the detection 
limit of that epoch. We use this selection method to be 
sensitive to short-term and long-term variability, as well 
as transients. This 5cr selection is quite conservative, and 
requires variability amplitudes increasing from |Ato| > 
0.1 mag for (m) ~ 18 mag up to |Ato| > 1.0 mag for 
(to) ^ 23 mag. 

We make the following cuts to the 5a variable source 
sample to remove artifacts: 

i) We remove sources with pipeline artifact flags 
indicating window bevel reflections or ghosts 
from the dichroic beam splitter. 

ii) We remove the brightest objects, with (to) < 
18.0, due to the large area subtended by the 
PSF which causes uncertainty in the back- 
ground subtraction. 

iii) We select sources within a radius of 0.55 deg 
of the center of the field, in order to avoid 
glints and PSF distortions, which are more 
prominent on the edges of the image, from 
un-corrected spatial distortions of photons 
recorded by the detectors. 

iv) We do not include objects that are within 
1.5 arcmin of a to < 17 mag source, to avoid 
regions affected by the bright source's PSF 
and ghost artifacts. Ghost artifacts can ap- 
pear within 30 — 60 arcsec above and below 
a bright source in the Y detector direction. 
Ghosts are point-like, and thus can only be 
identified from their Y detector position rel- 
ative to a bright source. While ghosts do 
not usually appear in the GALEX pipeline 
catalogs, we apply this cut since our forced 
aperture photometry could mistake ghosts for 
transient sources. 

v) We veto objects for which in the epoch of 
maximum \mk — {m)\/a{{m), k) or maximum 
flux, the aperture flux ratio of the object has 
h/ h.s > R90, where /g is the 6.0 arcsec 
radius aperture flux, /s.g is the 3.8 arcsec 
radius aperture flux, and Rgo is the max- 
imum cumulative aperture flux ratio mea- 
sured for 90% of the sources in the reference 
source sample used to calculate a{{m),k) in 
that epoch. This cut removes fluctuations in 
the background due to reflections from bright 
stars just outside the field-of-view, as well 
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Fig. 5. — Left: Gallery of good and vetoed variable sources dur- 
ing their epochs of minimum ("low") and maximum ("high") flux. 
Green circle shows a 6 arcsec (4 pixel) aperture radius. Sources 
vetoed for fa/fs.s > -R90 shown in panels a and b, a source vetoed 
as a likely ghost within 1.5 arcmin of a m < 17 mag source shown 
in panel c, and a manually vetoed source shown in panel d. For 
each source the grayscale is linear, and is scaled to the peak of the 
source in its "high" state. 

as epochs where the PSF is distorted due to 
a degradation in resolution which sometimes 
occurs in the Y detector direction. This also 
vetoes cases when the pipeline shreds a source 
into multiple sources, and the source is de- 
tected as variable because the center of the 
aperture is off-center from the peak source 
flux. 

vi) Finally, we visually inspect all of the remain- 
ing variable sources to remove any remaining 
artifacts that passed through the cuts above. 

Figure [S] shows a gallery of good sources and vetoed 
variable sources from our automated cuts (for a bright 
reflection in panel a, an off-center source in panel b, and 
a likely ghost in panel c) and manual cuts (for a diffuse 
reflection in panel d) . Our final GALEX TDS 5a variable 
sample after the cuts listed above has a total of 1078 
sources. 

3.2. Variability Statistics 

We characterize the variability of each 5a variable UV 
source using several statistica l measures. We measure t he 
structure function (following Idi Clemente et al.l ()1996l )). 



ViAt) 



-(|Ato 



At 



(2) 



where brackets denote averages for all pairs of points on 
the light curve of an individual source with i < j and tj — 
ti = At. The 2 day cadence of the observations combined 
with the seasonal visibility of the fields results in a dis- 
tribution of time intervals between observations (shown 
in Figure in]) that fall into 6 characteristic timescale bins; 
A<2d = 2±0.5 d, Ai4d = 4±0.5 d, Atgd = 6±0.5 d, Aigd 
= 8±0.5 d, Atiy,. = 0.96±0.14 yr, and Atjyr = 1.96±0.04 
yr. We measure the structure function in these 6 bins. 
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Fig. 6. — Histogram of the time intervals between all pairs of ob- 
servations for the 42 fields in GALEX TDS. Hatched regions show 
the time intervals over which the structure function is measured 
for all of the sources. 

and define Sd to be the maximum value of the struc- 
ture function evaluated for At2d, Af4d, Atgd, and Atsd, 
and Syi to be the maximum value of the structure func- 
tion evaluated for A^iyr and At2yr- We also measur e 
the intrinsic variability as defined hy Ses ar et aTl ()2007l ). 

and the maximum ampli- 



cr,„f = _^2^ where 

and e = ^ELi^((™).fc)' 
tude of variability, max(|ATO|). 



3.3. UV Light Curve 

In order to flag possible transient UV events that may 
be associated with a SN or TDE, we differentiate be- 
tween stochastic variability and flaring variability. We 
identify flaring UV variability as sources that show a 
constant flux > 10 days before the peak of the light 
curve, and do not fade more than 2a below the faintest 
pre-peak magnitude (measured > 10 days before the 
peak) . This selection criteria is tailored to the NUV rise 
times observed in SNe jG czari ct al. 2008b; B rown et al.l 
2009[ iGezari et all 120101: iMilne et al . 2010; Ofck et alj 
201C ) and TDE candidates dGezariet al. 2006, 2008^ 
2009[) . We define constant pre-peak flux as a light 
curve with a reduced x^/ < 3 (x^ = X]fc=i(™fc ~ 
(m))'^ /{a{{m),k)'^ /{n — 1), where n is the number of 
epochs > 10 days before the peak). For those sources for 
which there are only upper limits > 10 days before the 
peak, xl is set to 1. Flaring sources with no detections 
before the peak are further labeled as transients. Sources 
with xi/ > 3, or that fade below 2a of the faintest mag- 
nitude measured > 10 days before the peak are labeled 
as stochastically variable. We flag 116 flares, 145 tran- 
sients, 595 stochastically variable sources, and remain 
with 222 sources with neither light-curve classification 
flag. In Figure [7] we show example light curves of sources 
flagged as stochastically variable ('V'), flares ('F'), and 
transients ('T'). 

4. HOST PROPERTIES 

4.1. Archival Optical Imaging Catalogs 

We first characterize the host properties of the UV vari- 
able sources using archival optical u, g, r, j, and z pho- 
tometry and morphology fro m matches to the SD SS Pho- 
tometric Catalog, Release 8 (jAihara et al.ll201ll) (mum 
22 mag), the CFHTLS Deep Fields Dl, D2, and D3 
(miiiii 26.5 mag) and Wide Fields Wl, W3, and 
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Fig. 7. — Example light curves of sources flagged as stochastically 
variable ('V'), flares ('F'), and transient flares ('T'). Epochs > 10 
days before the epoch of peak flux are circled in red, and the epoch 
of peak fiux is circled in blue. Thick dashed line indicates the mean 
flux > 10 days before the peak. 

W4 (miiin ^ 25 mag) merged catalogs version T0005 ^, 
and the SWIRE ELAIS Nl and CDF S Region catalogs 
(miini ~ 24 mag) (jSurace et al.l 1200^ . For sources with 
matches in multiple catalogs, we use the match from the 
deepest catalog. We convert the CFHTLS magnitudes to 
the SP SS system using the conversions in Reg nault et al.l 
(|2009D . and the SWIRE Vega magnitudes to the SDSS 
system using the transformations measured for stellar ob- 
jects available at the INT WFS web sit e We then 
correc t for Galactic extinction using the iSchlegel et al.l 
(11998') dust map values for E{B - V) listed in Table H 
Figure [D shows the overlap of the GALEX TDS fields 
with the available archival optical catalogs. 

4.2. Pan-STARRSl Medium Deep Survey 

The GALEX TDS fields overlap with the PSl 
MDS fields MDOl (PSJCMMLSS), MD02 (PS.CDFS), 
MD04 (PS_COSMOS), MD07 (PS_GROTH), MD08 
(PS.ELAISN1), and MD09 (PS_VVDS22H). The 
Pan-STARRSl observations are obtained through a set 
of five broadband filters, (gpi, rpi, ipi, zpi, and ypi). 
F urther information o n the passband shapes is described 
in lStubbs et"all (|2010H . The PSl MD fields are observed 
with a typical cadence in a given filter of 3 days, with 
an observation in the gpi and rpi bands on night one, 
in the ipi band on night two, and the zpi band on 
night three, with j/pi-band observations during each of 
three nights on either side of the Full Moon. Image 
differencing is performed on the nightly stacked images, 
reaching a typical 5a detection limit of ~ 23.3 mag 
per epoch in the gpi, rpi, ipi bands and ~ 21.7 mag 
in the ypi band. Image difference de tections f rom the 
PSl Image Processing Pipeline (IPP; iMagnieil ([2006)) 
and an independ ent pipeline hosted by Harvard/CfA 
(|Rest et al.l l2005l ) are internally distributed to the 
PSl Science Consortium as transient alerts for visual 
inspection and classification. 

^ http:/ /terapix.iap.fr/rubrique.php?id_rubrique=252 

1 www .ast. cam. ac.uk/'^wfcsur /technical/phot om/colours 
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Fig. 8. — Comparison of the star/galaxy separator used in the 
PSl MD catalogs (thick dashed Unes) to matches in the SDSS (left) 
and CFHT (right) archival catalogs in the PS.GROTH field as a 
function of i-band magnitude. 

Deep stacks of the multi-epoch observations were 
generated to provide deep imaging with a 5a point- 
source limiting magnitude of 24.9, 24.7, 24.7, 24.3, 
23.2 mag in the ^pi, rpi, ipi, zpi, and ypi bands, 
respectively, and typical seeing (PSF FWHM) of ~ 
1.4, 1.3, 1.0, 1.0, 1.0 arcsec in the 5 bands respectively. 
The magnitudes are in the "natural" Pan-STARRSl sys- 
tem, m = — 2.51og(flux)-|-m', with a relative zeropoint 
adjust ment m' made in ea ch band for each individual 
epoch (jSchlafiv et al. 1 120121) before stacking to conform 
to the absolute flux calib ration in the AB magnitude 
system (jTonrv et al.ll2012l ). We convert the PSl mag- 
nitudes to the SDSS system using the bandpass trans- 
forma tions measured for stellar SEDs in Tonrv et al.l 
(I2012D. and c o rrect for Galactic extinction using the 
ISchlegel et all (|1998[ ) dust map values for E{B — V) 
listed in Table [H We obtain morphology information 
from the PSl IPP output parameters in the ipi filter 
for the PSF magnitude (PSF_INSTJ4AG), the aperture 
magnitude (PSF_AP_MAG), and the PSF- weighted frac- 
tion of unmasked pixels PSF_QF, to define a point source 
or extended source as: IF PSF_INST_MAG-PSF_AP_MAG 
< 0.04 mag AND PSF_QF > 0.85 THEN class = pt IF 
PSF_INST_MAG-PSF_APJ4AG > 0.04 mag AND PSF_QF > 
0.85 THEN class = ext We calibrated these parameter cuts 
by comparing sources detected in both the PSl MDS 
and archival optical catalogs. Figure [8] shows the PSl 
star/galaxy separation criteria for 110,804 sources de- 
tected in both PSl and SDSS catalogs in the PS.GROTH 
field, and for 169,461 sources detected in both PSl and 
CFHT catalog in the PS.GROTH field, with i < 22 mag, 
the faintest magnitude for 96% of the optical hosts of the 
GALEX TDS sources, and the magnitude limit where 
all three catalogs are complete. Even though the CFHT 
catalogs are deeper than SDSS, they do not attempt to 
separate stars and galaxies for i ^ 21 mag, and classify 
all sources fainter than this magnitude as point sources. 
However, it is clear from both comparison plots, that the 
PSl criterion of PSF_INSTJ4AG-PSF_AP_MAG < 0.04 mag 
does an even better job of separating the locus of stars 
from galaxies than both catalogs down to i 22 mag. 

4.3. Archival Redshift and X-ray Catalogs 

We also take advantage of the many archival X- 
ray and spectroscopic catalogs available from the over- 
lap of the GALEX TDS survey with legacy survey 
fields. In the PS_CDFS field we us X-ray catalogs 
from the 0.3 deg^ Chandra Extended CDFS survey 
(iGiacconi et al.l l2002t iLehmer et al] 120051: IVirani et al] 
12006D, and redshift catalogs from the VIMOS VLT Deep 



Survey (VVDS) ()Le Fevre et al.l 120041) . and a compi- 
lation of redshift catalogs from GOODS and SWIRE 
In the PSJCMMLSS field we us e X-ray catalogs 
from the 5.5 deg^ XMM -LSS survey (jChiappetti et al.l 
120051 : iPierre et all l2007|), an d redshift catalogs from 
VVDS (|Le Fevre et al.l l2005h . In the PS.COSMOS 
field, we use X-ray catalogs from the 1.9 deg^ X MM- 
Newton Wide-Field Survey ()Hasinger et al.l 120071) and 
the 0.9 deg2 Chandra COSMOS survey (|Elvis et al.l 
[2009f) . and redshifts from t he Magellan COSMOS AGN 
survey (iTrump et al.l [20071 [20091 th e VLT zCOSMOS 
bright catalog (|Lillv et al.ll2007n2009D . and the Chandra 
COSMOS Survey catalog (ICivano et al.llMl) . In the 
PS.GROTH field we use X-ray cat alogs from the 0.67 
degfChandra Extended Groth Strip (" Nandra et"al]|2005l: 
iLaird et al. 2009) and reds hift catalogs froin the DEEP2 
Galaxy Redshift Survey (jNewman et al.l I2012D . For 
PS.ELAISN1 we use the X-ray catalog fro m the 0.08 deg^ 
Chan dra ELAIS-Nl deep X-ray survey ([Manners et al.l 
12003( 1. Figure [D shows the overlap of the GALEX 
TDS fields with the archival X-ray surveys. Finally, 
we also match the sources with the ROSAT All-Sky 
Bright Source and All- S ky Survey Faint Source catalogs 
(|Voges et al.l[T999l[2000l) . 

5. CLASSIFICATION 

We classify the GALEX TDS sources using a combi- 
nation of optical host photometry and morphology, UV 
variability statistics, and matches with archival X-ray 
and redshift catalogs. Table [3] summarizes the sequence 
of steps we use to classify the sources, which we describe 
in detail below. 

5.1. Cross-Match with Optical Catalogs 

We first cross-matched our 1078 GALEX TDS sou rces 
with the archival u,g,r,i,z catalogs described in H4.ll 
with a matching radius of 3 arcsec. This radius is rec- 
ommended for matches between GALEX and ground- 
based optical catalogs (Budavari & Szalay 2008), and 
corresponds to a spurious match rate of only 1 — 2% 
at the high Galactic latitudes of the GALEX TDS fields 
(|Bianchi et al.|[20lll) . However, we found that there was 
a population of "orphans" (no optical match within 3 
arcsec) that were detected in their NUV low-state, and 
had a match between 3 — 4 arcsec with an optically iden- 
tified quasar. Given the strong likelihood that these are 
real matches, we increased our matching radius to 4 arc- 
sec. We use the star/galaxy classifications from the cat- 
alogs to label sources as point sources (pt) or extended 
sources (ext). This results in 878/1078 optical matches 
(81%), with the majority of sources without matches in 
PS_CDFS, which is only partially covered by the archival 
catalogs. 

We then match the sources that do not have archival 
optical matches to the PSl MDS catalog described in l4.2[ 
this increases the number of sources with optical pho- 
tometry and morphology (albeit without the u band) to 
1057/1078 (98%). Figure [S] shows a histogram of the r- 
band magnitude of the optical hosts, and the NUV — r 
colors of the GALEX TDS sources in their low-state. 
The optical hosts have a distribution that peaks at r ^ 21 
mag, over 3 magnitudes brighter than the detection limit 

http;/ /www. eso.org/~arcttura/CDFS_master/index. html 
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Fig. 9. — Left: Histogram of r magnitudes of optical matches. 
Right: Histogram of NUV — r colors of optical matches for those 
sources detected in the NUV during their low-state (solid lines), 
and those sources with upper limits during their low-state (dashed 
lines). 

of PSl MDS, and NUV - r - 1 mag. Sources not 
detected in their NUV low-state are shown as upper 
hmits in the NUV — r color histogram, and peak at 
NUV - r > 2 mag. 

5.2. Orphans 

We visually inspected the PSl stack images at the lo- 
cations of the 21 sources with no optical matches, and 
confirm that they are true orphan events. Furthermore, 
all of the orphans are undetected in their low-state in the 
NUV, with upper limits of NUV > (22.3 - 23.1) mag. 
Thus the orphan hosts are likely distant stars or faint 
galaxies (i.e., dwarf galaxies or high-redshift galaxies) 
that are undetected during their low-state in the opti- 
cal and NUV. 

5.3. Color and Morphology Cuts 

We first use the color and morphology of the opti- 
cal hosts to classify the GALEX TDS 5a UV variable 
sources. We define quasars as sources with optical point- 
source hosts with 



u- g < 0.7 
-0.1 < .g - r < 1.0 



(3) 



in order to avoid th e stellar locus and white dwarfs 
(jRichards et al.l I2002D . Note that this color selection 
can be contaminated by catacylismic variable stars 
(CVs), which overlap in color-color space with the 
quasar sample. Indeed, two of the sources classified by 
color as quasars are in fact spectroscopically confirmed 
CVs (VVDS22H_MOS05-05 and ELAISN1.MOS15-02). 
VVDS22H.MOS05-05 is R0TSE3 J221519.8-003257.2, a 
confirmed cataclysmic variables with a dwarf-nova type 
spectrum. We observed ELAISN1_MOS15-02 with the 
APO 3.5m telescope Dual Imaging Spectrograph (DIS) 
on 2011 May 3 and detected broad Balmer emission lines 
from a Galactic source, characteristic of a CV/dwarf nova 
spectrum. While these sources stood out easily because 
of their extreme magnitude of variability of \Am\ > 4 
mag, shown in Figure [TOl there may be lower amplitude 
CV events that are still hiding in our quasar sample. 
However, given the lo w surface density o f CVs relative 
to quasars in the sky (jSzkodv et al.ll20lll ). the expected 
contamination rate is consistent with the two CVs iden- 
tified. 

We define RR Lyrae stars as sources with optical point 
source hosts with 
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Fig. 10. — Maximum NUV variability amplitude as a function of 
low-state NUV magnitude for sources classified as RR Lyrae (red), 
M dwarfs (yellow), quasars (blue) AGN (green), and CVs (cyan). 
Dashed line shows the median Scr error selection function used to 
select the variable sources. 



0.75 <u-g < 1.45 
-0.25 < 5 - r < 0.4 
-0.2 < r - i < 0.2 
-0.3<i-z<0.3 



(4) 



(|Sesar et al.ll2010D . Note that the color cuts for quasars 
and RR Lyrae require the u-band, which is not available 
for sources with PSl-only matches. However, we define 
M dwarf stars as point sources with 



-i > 0.42 
z > 0.24 
g < 22.2 
r < 22.2 
i < 21.3 



(5) 



(6) 

u-band data. 



(jWest et al.ll201lD . which does not require 
We classify stars on the main stellar locus as those with 



1.0 < It -5 < 2.25 
0.4 < 5 - r < 1.0 
-0.2 < r - i < 2.0 
-0.3 <i- z <1.0 



(7) 



modified from lYannv et all (|2009l ). This color and mor- 
phology selection results in 37 RR Lyrae, 53 M-dwarf 
flare stars, 17 stars, and 325 quasars. Figure [TT] shows 
the optical color-color diagram of the sources with optical 
point-source hosts and u-band data, and their classifica- 
tions as quasars, RR Lyrae, M dwarfs, and stars. 

5.4. UV Variability Cuts 

Figure [12] shows the structure function on timescales 
of days and years described in f}3] for the sources clas- 
sified as RR Lyrae and quasars in H5.'3[ In Figure [T2] 
we show the ratio of the NUV structure function on 
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Fig. 11. — Colors of archival optical matches to UV variable 
sources with point-like optical hosts (black points). Dashed blue 
line shows the region in color-color space used to define quasars 
from optical colors and morphology alone. Sources with classifica- 
tion are color coded as quasars in blue, RR Lyrae in red, and M 
dwarf stars in yellow, and main stellar locus stars in cyan. Sources 
with archival X-ray matches are circled in purple. 

timescales of years to days (S'yr/S'd)- While quasars 
demonstrate a wide range of S'yr/S'd, all RR Lyrae have 
Syr/Sd < 3. We use this UV variability property to relax 
our color constraints, and increase our photometric sam- 
ple of quasars to all sources with optical point-source 
hosts with Syr/Sd > 3. This is equivalent to a struc- 
ture function po wer-law exponent c u t of 7 > 0.2, where 
S(At ) PC Ar (lHooket al.1 fl99l iVanden Berk et al.l 
[2OOI ISchmidt et al.ll2010ir This structure-function ra- 
tio selection results in the classification of another 30 
quasars. Two additional sources with optical point- 
source hosts have archival quasar spectra, resulting in a 
final quasar sample of 358 . We define active galactic nu- 
clei (AGN), as sources with optically extended hosts that 
show stochastic UV variability (see have an X-ray 

catalog match, and/or an archival spectroscopic classifi- 
cation. Figure [14] shows the optical color-color diagram 
of the sources with optically extended hosts, and those 
classified as AGN. This results in a sample of 305 AGN. 
We also add archival spectroscopic classifications for 6 
stars. This yields a total of 776/1078 (72%) sources clas- 
sified as an active galaxy (quasar or AGN) or variable 
star. 

5.5. X-ray Sources 

The archival X-ray catalogs overlap with ~ 8.45 deg^ of 
the GALEX TDS survey area. Within this area, 81/89 
quasars, 92/105 AGN, and 8/9 M-dwarf stars are de- 
tected in the X-rays. In addition, there are 9 optical 
point sources with X-ray matches that are likely quasars 
and M dwarfs just outside the quasar and M-dwarf color- 
color selection regions. UV variability selection appears 
to be selecting a similar population of active galaxies 
and M dwarfs as X-ray detection, since ^ 90% of the UV 
variability-selected active galaxy and M dwarf sample is 
also detected in the X-rays. However, only 2% of all the 
X-ray sources (the majority of which are active galaxies) 
are detected as UV variable at the selection threshold of 
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Fig. 12. — Histogram of the NUV structure function of sources 
classified as RR Lyrae (red) and quasars (blue) on a timescale of 
days and years. Red arrow and Blue arrow indicate the mean 
of Syr measured in the SPSS r-band for RR Lyrae and quasars, 
respectively from ISchmidt et al.l II2010I ). 
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Fig. 13. — Histogram of the ratio of the NUV structure function 
on timescales of years and days for sources classified as RR Lyrae 
(red) and quasars (blue). Dotted line shows the selection criteria of 
Sy/S^ > 3 used to classify sources with optical point source hosts 
as quasars. 

the GALEX TDS catalog. 

5.6. GUVV Catalog 

We also cross-match our GALEX TDS sample with the 
first and second GALEX Ultra violet Variability C ata- 
logs ( GUVV-1 fc GUV V-2) from I Welsh eTall (l2005l ) and 
Whe atlev et al.l (|2008D . These catalogs include 894 UV 
variable sources (Am > 0.6 mag in the NUV) selected 
from an analysis of archival GALEX AIS, MIS, DIS, 
and Guest Investigator (GI) fields with repeated observa- 
tions. With a cross-matching radius of 4 arcsec, we find 
a match with 36 GUVV sources. For the 15 matches that 
have GUVV classifications, are all classified by GUVV as 
active galaxies (AGN or quasars), which are in agreement 
with our GALEX TDS classifications. Of the 21 matches 
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Fig. 14. — Colors of archival optical matches to UV variables 
with extended optical hosts. Sources classified as AGN (by either 
stochastic U V variability, archival spectra, and / or an X-ray match) 
are color coded in green. Sources with matches with archival X-ray 
matches are circled in purple. Note that all galaxy sources with an 
archival X-ray match are defined as AGNs. 

without GUVV classifications, we find 5 sources classi- 
fied by GALEX TDS as RR Lyrae, 9 classified as active 
galaxies (AGN or quasars), 6 with optical point-source 
hosts, and 1 with a galaxy host. 

5.7. Unclassified Sources 

The remaining 302 unclassified sources include 91 with 
optical point-source hosts, which are likely stars, quasars 
with non-standard colors (high-redshift or reddened), or 
unresolved galaxies, 190 with galaxy hosts, and 21 or- 
phans. The 190 galaxy hosts may either be faint AGN 
with poorly constrained UV light curves, or hosts of UV- 
bright extragalactic transients. In Figure [T5] we show the 
maximum |Am| in the NUV as a function of low-state 
NUV magnitude of the remaining unclassified sources. 
The unclassified UV source with the most extreme am- 
plitude, ELAISN1.MOS15-09 with |Am| > 4.2 mag, was 
spectroscopically confirmed to be from the nucleus of an 
inactive galaxy si z = 0.1696, and its UV/optical fiare 
detected by GALEX TDS and PSl MDS was attributed 
to the tida l disruption of a st ar around a supermassive 
black hole (|Gezari et al.ll20f^ . Also in this sample is a 
UV transient spectroscopically confirmed to be a Type 
IIP SN 2010aqat z=0.086 (COSMOS_MOS26-29), whose 
UV/optical light curve from GALEX TDS and PSl MDS 
was fitted with early emissio n following SN sho ck break- 
out in a red supergiant star (jGezari "etalllMol) . Both of 
these spectroscopically classified extragalactic transients 
are labled in Figure [T5l 

Our 5(7 selection criteria translates to a limiting sen- 
sitivity to transients in a host galaxy with a magnitude 
Whost of a magnitude of 



mtr 



-2.51og(10- 



10^ 



(8) 



which ranges from mtrans 20.0 mag for mhost = 18 mag 
to mtrans '~ 22.7 mag for rrihost = 23 mag. Thus, our vari- 
ability selection threshold is less sensitive to transients 
in host galaxies with bright NUV fiuxes. On the red se- 
quence of galaxies, where Mnuv ~ —14.5 ()Wvder et al.l 
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Fig. 15. — Maximum NUV variability amplitude for sources 
without a classification. Dashed blue line shows the median Scr 
variability selection function as a function of mean magnitude. The 
nature of the optical host is indicated (point source, galaxy, or 
orphan), and flaring sources (including transients) are marked in 
yellow. 



I2007f ). this selection effect is not as much of an issue, 
since already for z > 0.05 one gets mhost > 22 mag. 
However, star-forming galaxies on the blue sequence are 
2.5 mag brighter in the NUV, and thus the host galaxy 
brightness can be a factor in reducing the senstivity to 
faint transients. For example, our GALEX TDS 5a sam- 
ple does not include SN 2009kf, a luminous Type IIP SN 
in a star-forming galaxy at z = 0.182 w hich we reported 
our GALEX TDS detection of in lBotticella et al.l (|20ldh . 
This source varied at only the 4.25(J level in the NUV 
during its peak. However, because this transient was 
selected from a spatial and temporal coincidence with 
a PSl transient alert, we could lower our threshold for 
variability selection in the UV. The systematic selection 
of SN and TDE candidates from the joint GALEX TDS 
and PSl MDS transient detections will be presented in 
future papers. 

6. DISCUSSION 

6.1. Classification Demographics 

Figure shows a pie diagram of the source classifica- 
tions. Out of the total of 1078 GALEX TDS sources, 
62% are classified as actively accreting supermassive 
black holes (quasars or AGN), and 10% as variable and 
fiaring stars (including RR Lyrae, M dwarfs, and CVs). 
Note that the relative fraction of the different classes of 
sources is sensitive to both their intrinsic magnitude dis- 
tribution, and the magnitude-dependent variability selec- 
tion function of the sample. Table 2] gives the GALEX 
TDS catalog, sorted by decreasing NUV amplitude, with 
the GALEX ID, R.A., Dec, low-state NUV magnitude, 
maximum amplitude of NUV variability (lAmmaxI), in- 
trinsic variability ((Ti„t), the structure function on day 
(^d) and year (Sy) timescales, the characteristics of the 
NUV light curve: fiaring (F) or stochastically variable 
(V), the morphology of the matching optical host: point- 
source (pt) or extended (ext), the color classification of 
the matching optical host: RR Lyrae (RRL), M dwarf 
star (Mdw), star (Star), or quasar (QSO), the archival 
redshift, an X mark if there is a match with an archival 
X-ray source, and finally the GALEX TDS classifica- 
tion: RR Lyrae (RRL), M dwarf star (Mdw), star (star), 
quasar (QSO), AGN, UV fiaring source or UV transient 
source with galaxy host (Galaxy Flare or Galaxy Trans), 
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Fig. 16. — Pie chart of GALEX TDS classifications: sources clas- 
sified as stars (Star), M dwarfs (Mdw), RR Lyrae (RRL), quasars 
(QSO), and active galactic nuclei (AGN), and sources with no clas- 
sification with galaxy hosts (Galaxy), no host (Orphan), and point- 
source hosts (Point). 

UV flaring source or transient source with point-source 
optical host (Point Flare or Point Trans) or UV flar- 
ing source or transient source with orphan optical host 
(Orphan Flare or Orphan Trans), stochastically vari- 
able source with a point-source optical host (Point Var), 
stochastically variable orphan (Orphan Var), or none of 
the above (?). 

In Figure [17] we show the cunimulative surface density 
distribution of classifled UV variable sources as a func- 
tion of peak magnitude (high(NUV)) and maximum am- 
plitude {ma.x{\ANUV\)). For the variable UV sources, 
these correspond to total surface densities of 8.0 ± 3.1, 
7.7 ± 5.8, and 1.8 ± 1.0 deg"^ for quasars, AGNs, and 
RR Lyrae, respectively. For the transient source, we can 
calculate a total surface density rate, #/(area x tee), 
where t^s is the effective survey time at the cadence that 
matches the characteristic timescale of the transient. For 
extragalactic transients such as young SNe and TDEs, 
which vary on a timescale of days, we use the time in- 
tervals for which the fields were observed with a cadence 
of 2.0 ± 0.5 days. If we include all flaring or transient 
GALEX TDS sources with a galaxy host, this yields a 
surface density rate of 52 ±38 deg~^ yr~^ for extragalac- 
tic transients. For M dwarfs which vary on timescales 
shorter than an individual observation, we use the total 
exposure time for each epoch. If we assume a survey with 
a cadence of 2 days and iexp = 1-5 ksec, this translates 
to a surface density rate for M dwarfs of 15 ± 10 deg~^ 
yr-i. 

6.2. UV Variability Properties of Classified Sources 

Various optical studies of rest-frame UV variabil- 
ity in high redshift quasars have demonstrated that 
the variabil ity of AGN increases with decreasing rest 
wavelength ("di Clement e et aLlll996t IVanden Berk et all 
[2004; Wilhitc et al. 200l) In Figure [H we show his- 
tograms of (Tint for the UV variable sources with clas- 
sifications. Quasars show a distribution of dint with a 
mean that is 2 times larger than measured at opti- 
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Fig. 17. — Cummulative distribution of surface density of 
UV variable sources with GALEX TDS classifications; M dwarfs 
(Mdw), RR Lyrae (RRL), quasars (QSO), active galactic nuclei 
(AGN), and extragalactic transients (GAL-Flare) as a function 
of peak magnitude (high(NUV)) (left) and maximum amplitude 
(max(|AAri/y|)) [right). 

cal wavelen gths from the SDSS Stripe 82 sample from 
ISesar et al.. (2007). This effect is even more pronounced 
in the magnitude of the structure function on years 
timescales (Sy), which has a mean that is 5 times larger 
than the mean mea sured in the r-band (Ae// — 6231) 
from [Schmidt et all (2010). This trend is consistent with 
the wavelength dependent rise in variability amplitude 
observed in the structure function for quasars in the rest- 
frame UV ([Vanden Berk et al.[ [20040 and observed UV 
([Welsh et al.ll2011[ ). 

The fact that AGN become bluer d uring high states of 
fiux ( Giveon et al.[[1999t [Geha et al.l [20031: [Gezari et all 
2008 a|) has been attributed to increases in the charac- 
teristic temperature of the accretio n disk in response to 
increases in the mass accretion rate ([Perevra et al.| [2006[: 
[Li fc Cad [20081 ). However. [Schmidt et al.[ ([2012[ ) argue 
that the color variability observed in individual quasars 
in their SDSS Stripe 82 sample is stronger than expected 
from just varying the accretion rate (M) in accretion 
disk models. In a future study, we will use simultaneous 
UV and optical light curves from GALEX TDS and PSl 
MDS for our 358 individual quasars to test this result 
with a larger dynamic range in wavelength. 

For the subsample of 95 quasars with archival red- 
shifts (zmean — 1-26, (T^ — 0.39), in Figurc [19] wc plot 
o'int VS. the low-state NUV absolute magnitude, and 
find a steep negative correlation fitted by log(crint) — 
(1.6 ± 0.1) + ^Mnuv, where /S = 0.24 ± 0.04, in ex- 
cellent agreement with the trend for increased variabil- 
ity in lower luminosity quasars see n from optical ob- 
serva tions with /? = 0.246 ± 0.005 ([Vanden Berk et "all 
[20041 ) , and shallower tha n expected for a Poissonian pro- 
cess which has P = 0.5 ([Cid Fernandes et al.|[2000l ). We 
also show the subset of 68 AGN with archival redshifts 
(•Zmcan = 0.64, (T^ = 0.55), which clearly do not show a 
relation between tJint and low-state NUV absolute mag- 
nitude. This is most likely a result of dilution of the 
variability amplitude from the contribution of the host 
galaxy in the NUV. 

The largest values of |Am| (plotted in Figure [TU[) are 
found in RR Lyrae and M dwarfs, with a tail of large 
amplitude variations reaching up to |Ato| = 2.9 mag in 
RR Lyrae and up to | Am| — 4.6 mag in M dwarfs. In the 
optical, the RR Lyrae structure function is very weakly 
dep endent on times cale, with an amplitude of 0.1 — 0.2 
mag [Schmidt et all (WlO). The NUV structure function 
also shows a weak dependence of amplitude on timescale 
when comparing the structure function on days to years 
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Fig. 18. — Intrinsic NUV variability as a function of low-state 
NUV magnitude for sources classified as RR Lyrae (red), M dwarfs 
(yellow), quasars (blue) and AGN (green). Arrows shows the me- 
dian (Tint measured in the optical for quasars (blue arrow) , RR 
Lyrae (red arrow), and M dwarfs (orange arrow) from lSesar et aO 
||2007|) . 
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Fig. 19. — Intrinsic NUV variability as a function of low-state 
NUV absolute magnitude, Mjvc/v = low{NUV) - DM where DM 
is the distance modulus, for the subsample of quasars (blue dots) 
and AGN (green circles) with catalog redshifts. Dashed blue line 

shows the fit to the quasars to log(a-i„t) oc The spectro- 

scopically classified extragalactic transients (TDE PSl-lOjh and 
SN 2010aq) are labeled. 

timescales, however, with an amphtude that is ~3 times 
larger in the NUV than in the optical This wavelength 
dependence on variability amplitude can be explained if 
the variability is driven by variations in surface tempera- 
ture from pulsations of the stellar envelope, where higher 
states of flux are associated with higher temperatures 
(jSesar et al.ll2007D . 

7. CONCLUSIONS 

We provide a catalog of over a thousand UV variable 
sources and their classifications based on optical host 
properties, UV variability behavior, and cross-matches 



with archival X-ray and redshift catalogs. This yields a 
sample of 53 M dwarfs, 37 RR Lyraes, 358 quasars, and 
305 AGN. We find median intrinsic UV variability am- 
plitudes in RR Lyrae and quasars that are factors of > 3 
larger than at optical wavelengths, consistent with the 
expectation for higher temperatures during higher states 
of fiux. The regular cadence and wide area of GALEX 
TDS enables us to systematically discover persistent and 
transient (i.e. tidal disruption of a star) accreting super- 
massive black holes over wide fields of view, study the 
contemporaneous UV and optical variability of variable 
stars, and catch young core-collapse supernovae within 
the first days after explosion. The overlap of the GALEX 
TDS with PSl MDS and multiwavelength legacy survey 
fields will continue to be helpful for classifying transient 
sources in these heavily observed fields. We also measure 
the surface densities of variable sources and the surface 
density rates for transients as a function of class in the 
UV for the first time. 

With GALEX TDS we are only scratching the surface 
of UV variability. Our 5<j sample of 1078 sources is less 
than 0.3% of the 419,152 UV sources in the field (with 
an average density of 1.1 x 10"* UV sources per square 
degree down to miim = 23 mag). Looking to the future, 
the discovery rate for UV variable sources and UV tran- 
sients could increase by several orders of magnitude with 
the launch of a space-based UV mission with a wide field 
of view (several deg^), a survey strategy of daily cadence 
observations over ~ 100 deg^, and detectors with an or- 
der of magnitude improved photometric precision relative 
to GALEX. In the optical sky, 90% of quasa rs vary with 
Ojnt > 0.03 mag on the timescales of years (jSesar et al.l 
|2007| ). Given the factor of ~ 2 larger crint observed for 
quasars in the NUV, one could achieve a nearly complete 
sample of low-redshift quasars with photometric errors of 
a{m) 0.01 mag. In coordination with a grou nd-ba sed 
optical survey, such as Pan-STARRS2 (|Burgeti[20ll l or 
LSST this could yield the simultaneous UV and op- 
tical detection of « lO'* variable quasars and « lO'^ RR 
Lyrae and M dwarfs, as well as increase the discovery 
rate of UV-bright extragalactic transients (young SNe 
and TDEs) by a factor of « 100. 
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of Hawau, the Pan-STARRS Project Office, the Max- 
Planck Society and its participating institutes, the Max 
Planck Institute for Astronomy, Heidelberg and the Max 
Planck Institute for Extraterrestrial Physics, Garching, 
The Johns Hopkins University, Durham University, the 
University of Edinburgh, Queen's University Belfast, the 
Harvard-Smithsonian Center for Astrophysics, and the 
Las Cumbres Observatory Global Telescope Network, In- 
corporated, the National Central University of Taiwan, 
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GALEX Time Domain Survey 



TABLE 1 
GALEX TDS and PSl MDS 



Survey 


Field of View 


Plate Xcale 


PSF FWHM 


™lim per Epoch 


Cadence 


Seasonal Visibility 




(deg) 


(arcsec/pixel) 


(arcsec) 


(mag) 


(days) 


(months) 


GALEX TDS 


1.1 


1.5 


5.3 


23.3 


2 


~ 1 


PSl MDS 


3.5 


0.258 


1.0 


23.0 


3 


~ 6 



TABLE 2 
GALEX TDS Fields 



Name R.A. (J2000) Dec (J2000) E{B - V) Epochs 

(deg) (deg) (mag) 



o _A IVi iVl o o _iVi U o U U 


OO 


.OoU 


Q 

— O. 


1 /I n 
, 14U 


U 


noQ 
.UZo 


ov 
Z / 


o _A1\ 1 iVi LlOO _iVi U o U i 




.oUU 


Q 

— O. 


/inn 


U 


.UZo 


Z f 


"DC VA/rA/TT GC AJ'OG/IO 


35. 


.000 


—3. 


.950 


0. 


.020 


14 


T3C3 YTV/TA/TT OC3 TV/Ti^GflQ 

r o -AiVlIVl J_(0 o -iVLVJ o Uo 


OO. 


.tsio 


A 

— 4. 


ocn 


U. 


.UZo 


07 
Z 1 


Jr o _AiVLIVl J_iO o _iVLU o U4 


"in 

OD. 


.yuu 


A 

— 4. 


.4ZU 


U. 


.UZo 


07 

Z r 


T3G YTV/TA^T GG A/r/~\Gnc; 


OO. 


.zUU 


c 

— o. 


.UoU 


U. 


noo 
.UZZ 


Ofi 

ZD 


T3G Y1\/T"\TT GG A/r/~\Gnf; 
r o -AiVilVl J_iO o -iViiJ £) Uu 


oD. 




c 

— o. 


onn 
.ZUU 


U. 


no? 
.uz / 


O/t 

Z4 


T3G OTM7G A/Ti^Gnn 


CO 
OO. 


.100 


—27. 


.800 


0. 


.008 


11/1 
114 


T3G A/T/^Gm 


oz. 


.Uiz 


—zo. 


oi o 
.ZIZ 


U. 


.UUo 


Qn 
oU 


"DG OTMT'G A/rr~iGno 


CQ 

OO. 


, iz4 


— ZD. 


Qno 


U 


.uuy 


OO 

zy 


PG r'FtTT'G A/TMGri'^ 


04, 


1 fit; 
. iuo 


— Z / . 


QT O 
,OlZ 


U, 


m o 
,U1Z 


oU 


r>G (^"TlTTG A/TOGri/l 


OZ, 


m n 

.yiu 


OQ 

— Zo. 


.oUU 


U, 


.uuy 


Qn 
oU 


r>G (^"TlTTG A/TOGnc; 


oz, 


111 
,111 


— z / . 


.Z / D 


U, 


m n 
.UlU 


oU 


r>G (^TlTTG A/f OGHfi 


OO 


.y /U 


OQ 

— Zo. 


.004 


U, 


m n 
.UlU 


D 


o J o 1\ i o _iVi \JoZi. 


loU. 


.oUU 


+0, 


1 nn 
. lUU 


U, 


noQ 
.UZo 


1 C 
10 


r 0_0'vyolviV_/0_iVi WoZZ 


149 


.OUU 


-f-o. 


1 no 

. lUU 


n 

u. 


097 
.uz i 


1 


PS.COSMOS.MOS23 


151, 


.000 


+2, 


,200 


0, 


,024 


24 


PS.COSMOS.MOS24 


150, 


.000 


+2, 


,200 


0, 


,020 


26 


PS.COSMOS.MOS25 


149, 


.000 


+2, 


,300 


0, 


,023 


13 


PS.COSMOS.MOS26 


150. 


.500 


+1, 


,300 


0, 


,023 


26 


PS.COSMOS.MOS27 


149. 


.500 


+1. 


.300 


0, 


.019 


27 


PS.GROTH.MOS01 


215. 


.600 


+54, 


.270 


0, 


.011 


17 


PS.GROTH.MOS02 


213. 


.780 


+54, 


.350 


0, 


.015 


16 


PS.GROTH.MOS03 


214. 


.146 


+53. 


.417 


0, 


.009 


17 


PS.GROTH.MOS04 


212. 


.400 


+53. 


.700 


0, 


.011 


18 


PS.GROTH.MOS05 


215. 


.500 


+52, 


.770 


0, 


.008 


19 


PS.GROTH.MOS06 


214. 


.300 


+52, 


.550 


0, 


.008 


8 


PS.GROTH.MOS07 


212. 


.630 


+52, 


.750 


0, 


.009 


19 


PS.ELAISN1.MOS10 


242. 


.510 


+55, 


.980 


0, 


.007 


17 


PS.ELAISN1.MOS11 


244. 


.570 


+55, 


.180 


0, 


.009 


17 


PS.ELAISN1.MOS12 


242. 


.900 


+55, 


.000 


0, 


.008 


18 


PS.ELAISN1.MOS13 


241, 


.300 


+55, 


.350 


0, 


.007 


18 


PS.ELAISN1.MOS14 


243. 


.960 


+54, 


.200 


0, 


.010 


19 


PS.ELAISN1.MOS15 


242. 


.400 


+54, 


.000 


0, 


.011 


21 


PS.ELAISN1.MOS16 


241. 


.380 


+54, 


.450 


0, 


.010 


19 


PS.VVDS22H.MOS00 


333, 


.700 


+1, 


.250 


0, 


.040 


39 


PS.VVDS22H.MOS01 


332. 


.700 


+0, 


.700 


0, 


.046 


35 


PS.VVDS22H.MOS02 


334. 


.428 


+0, 


.670 


0, 


.057 


38 


PS.VVDS22H.MOS03 


333. 


.600 


+0, 


.200 


0, 


.058 


27 


PS.VVDS22H.MOS04 


334. 


.610 


-0, 


.040 


0, 


.093 


24 


PS.WDS22H.MOS05 


333. 


.900 


-0, 


.720 


0, 


.102 


35 


PS.WDS22H.MOS06 


332. 


.900 


-0, 


.400 


0, 


.113 


33 



14 



Gezari et al. 



TABLE 3 
GALEX TDS Classifications 



Archive PSl Classification 



Step 


N 1 

^ ' unclass 
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ext 


pt 


ext 


orphan 
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RRL 
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Star 


AGN 
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QSO color cut 
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RRL color cut 
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Mdw color cut 
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stellar locus color cut 
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Syr/ > 3 
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X-ray /spec match 
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37 








302 










21 


358 


37 
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22 
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93 
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TABLE 4 
GALEX TDS Catalog^ 











NUV 












Optical 




X-ray 


Class 


ID 


R.A. 


Dec 


nijciTO 


|Ammax| 


O'int 


Sd 




LC 


Morph 


VAB 


Color z 






GROTH.MOSOl-21 


216.1622 


54.0911 


22.54 


4.60 


1.04 


0.80 


1.13 


V 


pt 


14.92 


Mdw 




Mdw 


VVDS22H.MOS05-05 


333.8326 


-0.5491 


21.14 


4.47 


0.99 


0.90 


0.74 


F 


pt 


21.24 


QSO 




CV 
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F 
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* Table 4 is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its format and content. 



